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I. INTRODUCTION 



A. MOTIVATION OF STUDY 

Microwave search radars installed on ships, or at shore stations, operating close 
to the sea surface, have a strong interaction with the propagating medium and the sea. 
The reflected field is a function of the transmitting frequency, the incidence angle, and 
the electric properties of the surface boundary. The atmospheric medium also interferes 
with the electromagnetic wave by refraction, absorption, and attenuation. 

The present study tries to evaluate quantitatively the effects of sea roughness 
and special atmospheric conditions, like ducting, on electromagnetic propagation. It is 
shown that the propagation is considerably influenced by the media between the radar 
and the target. 

The motivation of analyzing the irregularities of the sea surface with periodic or 
random variations, together with atmospheric refractivity, is to predict and optimize 
the detection ranges of a search radar. 

B. SCOPE OF WORK 

This study begins by introducing the Rayleigh criterion and the Fresnel zone 
model for discriminating between a rough and a smooth surface. The specularly re- 
flected wave interferes with the direct propagating wave, resulting in a constructive 
and destructive interference. 

Computation and illustration of detection patterns, depending on radar parame- 
ters and sea surface characteristics, are performed. The resulting plots show a strong 
effect of constructive and destructive interference between the direct and reflected fields. 
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It is shown that under some special conditions the detection range could be doubled 
relatively to free space range. 

The abnormal propagation caused by strong vertical gradient of tropospheric 
refractivity, results in a considerably extended horizontal detection range. From con- 
servation of energy point of view, the increased density of energy through one path 
results in a decreased density through other paths causing “holes” in detection. 

The results of this study show that by appropriate design and operation, detection 
ranges of a radar could probably be improved. Also, the limited detection range could 
be predicted and optimized. 

C. ORDER OF WORK 

Chapter II presents a theoretical study on surface roughness, specular reflection 
and diffused scattering. The Rayleigh criterion and Fresnel zones for surface roughness 
are introduced. 

Chapter III introduces a simplified model for propagation above a spherical sur- 
face. The divergence factor and the reflection coefficient from the sea surface axe 
discussed in detail. The shadowing effect for radar with very low grazing angles, and 
the depolarization by the reflecting surface are also introduced. 

Chapter IV is based on the study of Chapter II and III. This chapter discusses 
the lobing phenomena. Computer plots of detection contours are illustrated. Conclu- 
sions regarding improved detection range for low flying aircraft and surface targets axe 
suggested. 

Chapter V introduces abnormal propagation in the atmosphere. The atmospheric 
conditions for this phenomena to occur, and the effect on radar detection ranges axe 
discussed. 



2 



Chapter VI is the summary and conclusion of this study. The Appendix gives an 
analytical derivation of Fresnel zones. 
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II. RAYLEIGH CRITERION AND FRESNEL 
ZONES FOR SURFACE ROUGHNESS 



A. INTRODUCTION 

When an electromagnetic short wave (VHF and above) propagates above the 
earth’s surface, it is affected by the earth’s surface and its atmosphere. The reflected 
field of the electromagnetic wave from such boundaries depends on the angle of in- 
cidence, the wavelength and the characteristics of the reflecting plane (conductivity, 
permittivity, and permeability). The difficulty of qualitatively analyzing the influence 
of the medium on the propagating wave is due to the fluctuating irregularities of the 
surface boundaries. These irregularities impose limitations on the performance of com- 
munication and radar systems. These random fluctuations of the parameters, defining 
a practical terrain profile, are the main obstacle in establishing its relationship with 
the properties of the electromagnetic wave. For this reason the prediction of radar 
and communication performance is limited. In other words, the effects on propagation 
introduces degradations in radar performance. 

In this chapter we will discuss the theory of specular reflection and its effect 
upon field strength on the assumption that the reflection coefficient of a perfectly 
smooth earth are known. Effects of surface roughness will also be discussed. 

B. ANALYTICAL MODEL AND CONSIDERATIONS 

1. Surface Roughness: General Definition 

The first stage in constructing a model for the propagation problem is to 
analyze and understand the reflecting boundary which is usually varying in position 
and time, and therefore is not plane. 
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The basic difference between a rough and a smooth surface is that a smooth 
surface will reflect the incident plane wave specularly in a single direction, whilst a 
rough surface will scatter the wave in various directions (see Fig 2.1, 2.2). From this 
definition one can conclude that the roughness of the surface is determined basically 
by the angle of incidence and the frequency of the plane wave. A smooth surface is a 
singular limiting case of a rough one. 



R 




Figure 2.1: Reflected Waves from Perfectly Smooth Surface [Ref. 1]. 

In Chapters III and IV we will introduce the sea medium which will include 
the time varying nature of the scattering surface and its effects such as the fading of 
the received signal by the shape and movement of the water waves [Ref. 3], 

2. Fundamental assumptions 

The empirical and analytical theory of rough surface scatter is usually sim- 
plified by some of the following assumptions [after Ref. 3]: 

1) Far field approximations only. 
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Figure 2.2: Reflected Waves from Rough Surface [Ref. 2]. 
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2) Shadowing effects are neglected. 

3) The radius of curvature of the scattering surface is much larger than the wave- 
length of the incident wave. 

4) Multiple scattering is neglected. 

5) Particular models of surface roughness are restricted (periodic and random). 

3. Roughness Criteria 

As mentioned earlier, the efficiency of specular reflection from a rough 
surface is a function of the ratio of the surface roughness dimensions to the wavelength, 
and of the grazing angle. We will try to define the values of the grazing angle, the 
wavelength, and the surface roughness at which the specular reflection changes into 
diffuse scattering. 

The best known criterion involving these parameters is that of Lord Rayleigh. 
The path difference between the two waves is 

8 = 2A/isin#. (2.1) 



The phase difference A (f> corresponding to this path difference is: 

. . , r 47 tA/i . . 

A<p = k8 = — - — sin 0, 



( 2 . 2 ) 



where A is the wavelength and the propagation number is given by k = y-- 

For small phase differences, there is almost no effect by the roughness and 
the surface behaves as a smooth surface. This can be achieved if ^ or 9 are very small 
which leads to A <f> « 0. Increasing ^ and/or 6 will increase the phase difference. 
When the phase difference is equal to 7 r the two waves will be in phase opposition 
and cancel in that direction. According to conservation of energy, the energy of the 
electromagnetic wave is scattered in other directions. Thus we conclude that for: 

1) A <f> = 0 the surface reflects specularly and is smooth . 

2) A(f> = 7r the surface scatters and is roug h [Ref. 1], 

A criterion for discriminating between rough and smooth surfaces is chosen 
arbitrarily as 7 t/ 2 which is the average phase difference between these two extremes 
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[Ref. 5]. Thus a surface is considered to be smooth for the height of roughness Ah to 
be smaller than: 



Ah < 



8sin0 



(2.3) 



This criterion, sometimes called the Rayleigh criterion, is a simplified argument based 
on ray theory for describing the irregularities of a surface. 

The Rayleigh model excludes, by definition (because of the ray assumption), 
the case when the roughness of the surface elements are comparable in size to the 
wavelength of the propagating wave. 

The sea surface in such a case usually consists of irregular distributions of 
sizes and shapes of the waves. An accurate solution in this case requires complex 
methods not yet developed to satisfy boundary conditions of a random rough surface. 
However, in the literature, simplifications have been made in order to calculate the 
effective reflection coefficient, using sea surface models with sinusoidal or wedge shaped 
waves [Ref. 3]. 

4. Active Scattering Region 

As mentioned above, the complexity of the reed surface is so severe that no 
analytical model can describe it accurately. In this section we will try to consider the 
effects of scattering in terms of a simple model called the Fresnel zones, or the glistening 
surface. Ray analysis suggests that reflections would be confined to a small region of 
the rough surface, which depends on the transmitter receiver configuration. 

The glistening surface, or the Fresnel near zones, is defined as that part of 
a rough surface that reflects a significant amount of electric field of the propagating 
wave into the receiver, at given heights of the source and the receiver [Ref. 4]. 

Applying Huygen’s principle to the ray theory, one should interpret the re- 
flected ray from a single point on the surface as the fact that an energy source illumi- 
nating a large portion of the surface creates currents which radiate in all directions. At 
any point in space the total field is the sum of the radiation from these currents and 
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Figure 2.3: Glistening Zones on a Reflecting Plane [Ref. 4]. 

the direct propagating field. Figure 2.3 illustrates a scenario consisting of a transmitter 
at height z\ and a receiver at height operating above a smooth plane. 

The direct ray path is Rd and the indirect path from a point on earth is 
R\ + i? 2 • As mentioned, the incident field induces currents on the plane xy which radiate 
a reflected field. We want to find the zones in which the path difference between the 
direct and the indirect path is constant and equals to A<f>. The path difference is [Ref. 
1,3]: 



R i R 2 — Rd — 8. 


(2.4) 


Since the direct path is unchanged and 8 is constant: 

Ri -f- R 2 = Const. = Rd 4" 8 


(2.5) 


Putting the path difference 8 into phase difference notation: 

A <{> = 8 • k, 


(2.6) 


O 

*1 

t> 

II 

^ to 


(2.7) 
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Equation 2.5 is an equation of ellipsoids with T and R as foci and the direct path 
Rd as axis of revolution (see Figure 2.3). The geometric locus of all points on the xy 
plane which intersect with the ellipsoids tire ellipses which describe a constant phase 
difference, A <f>. In order to describe more zones of other constant path difference we 
increase S, for instance, by A/2, or 7r radians (for phase interpretations). 

We find a family of ellipses which describes zones of constant phase. Each 
zone differs from its adjacent zone by 7 r, so the radiation from two successive zones 
is in opposite phase. They do not cancel totally because of the slow variation in the 
amplitude between zones. We shall see that the first zone, n = 1, is the contributing 
zone [Ref. 4]. 

a. Fresnel zones dimensions 

In order to apply smooth earth phase difference formulas one should 
know the dimensions of the Fresnel zones, which can give an indication as to what the 
minimum dimensions of a smooth area should be for these formulas to be valid. 

The approximate path difference between the direct path and the geo- 
metrical point (only one) of reflection is: 

So « — ; ( 2 . 8 ) 

r 

the nth ellipse path difference is determined by: 

_ „ n A .. 

S n & Sq + — . (2.9) 



By using this approximation and substituting it into Equations A. 13, 
A. 15, and A. 11 and using the assumptions of Equation A. 16 (Appendix A), we find 
useful formulae for the semiminor axis (yi„), the semimajor axis (xi n ) and the center 
(x 0 „) of the nth ellipse: 



Vln 




N 




1 , (^ 1 + Z ?) 2 ’ 

1 ^ nAr 



r 

Zln ~ XQn ± - 



\/i + S 

1 , (*l+*2) 2 

1 ' nAr 



A 

^1 n — ^0 n 



± 



( 2 . 10 ) 

( 2 . 11 ) 
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and 



ZOn ~ X 



rl + 2il ^j 



2 1 + (*> + ? 2 ) 2 . 

1 nAr 

The area of the nth ellipse is: 

A n = 7Tfnyin- 

By substituting Equations 2.10 and 2.11 into Equation 2.13, we obtain: 

7r r VnXr(l + & ) 

"" 4 (i + ^f 2 ) 3 ' 2 ' 



( 2 . 12 ) 



(2.13) 



(2.14) 



b. Special Practical Cases 

The minimum area for smooth surface assumption (Fresnel zone) calcu- 
lated by Equation 2.14 can be simplified for special practiced cases of a radar tracking 
an aircraft: 



4zi 

z-i 



« 1. 



(2.15) 



The height of the radar antenna is z\ and the height of an aircraft is 
z 2 . Since the first zone is the most contributing area to the total field received at R 
(Fig. 2.3), we should consider only this zone (n = 1). At the reflection point for the 
lowest lobe of the interference pattern, r ^ zi,z 2 and <5 = 



yn 



xu 



x oi 



lplified to [Ref. 4]: 


2x/2z! « *° , 

1 3^27, ’ 

A 


(2.16) 


.8 ^ 
Xoi± ^ , 


(2.17) 


12z? 


1 

A ‘ 


(2.18) 



Two cases are considered to illustrate the dimensions and positions of 
the first Fresnel zone. In the first case, the range, the height of the source z i, and the 
height of the receiver z 2 , are given and it is required to determine the position and 
dimensions of the first zone. In the second case r and z\ are given and we want to find 
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dimensions of the first zone. In the second case r and Z\ are given and we want to find 
the dimensions and position of the first zone which will produce the reflected wave to 
form the lowest lobe of the interference pattern for which 8 0 = j. The calculated values 
for given grazing angle and mean height of sea waves axe given in Table 2.1 [Ref. 3,4]. 



TABLE 2.1: FRESNEL ZONES AREA A 1 [REF. 4]. 



Parameter 


Case 1 


Case 2 


A [cm] 


10 


100 


10 


100 


r [m] 


10,000 


10,000 


10,000 


10,000 


z-i [m] 


1000 


1000 


2.5 


2.5 


z\ M 


100 


100 


100 


100 




20m 


20m 


5cm 


50cm 


Xoi [m] 


910 


943 


9350 


6830 


Vn M 


9.1 


29 


6.6 


44 


6 H 


83 


262 


615 


2760 


*[°] 


6.3° 


6.3° 


0.59° 


0.72° 


^max [m] 


0.11 


1.1 


1.22 


1.25 


A x [m 2 ] 


2,370 


23,900 


128,000 


382,000 



Table 2.1 indicates that the minimum areas of the first zone for smooth 
assumptions seem to be large, and are larger with increasing wavelength. For mean 
peak to peak wavelength of almost lm, in case 2, the smooth area for the first lobe of 
interference is much larger than the smooth area in case 1, especially for short wave 
lengths. 

The contribution of energy to the receiver from each zone is calculated 
by integrating the current excitation of each zone. Though the successive zones differ 
in phase by 7 r radians, it is not evident that their contribution to the received energy 
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cancels out totally and that only the first zone or part of it contributes to the received 
energy. 

The reason for this doubt could be caused by the variation in the am- 
plitude of excitation within a single zone, especially in small grazing angles. From later 
research it seems that the first Fresnel zone is likely to be the most important zone of 
contribution to the received energy, together with two other areas near the source and 
near the receiver which differ from the Fresnel zones. 

In this discussion only the phase of the reflected radiation has been con- 
sidered. We did not consider the conductivity and dielectric constant nor the reflection 
coefficient. These will be discussed in the next chapter. The sea surface is more com- 
plex than the land surface, with time variations in addition to fluctuations in amplitude 
and phase of the scatterers. 



13 



III. EARTH’S CHARACTERISTICS 



AFFECTING PROPAGATION 

A. INTRODUCTION 

After discussing the roughness of the surface in Chapter II and deriving roughness 
criteria and a model, we want to discuss in this chapter the major characteristics of the 
earth’s surface which will affect the magnitude, phase, and direction of a propagating 
wave. 

Firstly, we will consider the effects of a spherical surface and refraction effects and 
present a model for horizon propagation range, taking into account the convex shape 
of the earth’s surface. Secondly we will analyze the reflectivity of the surface. 

Compared with propagation in free space, the presence of terrain affects the Radar 
Cross Section (RCS) of scatterers. The reflection from terrain creates an interference 
pattern which will be discussed in detail in Chapter IV. Also, the terrain casts a shadow 
which introduces diffraction phenomena. The main outcome will be a reflection coef- 
ficient for a rough and a smooth surface, which will define quantitatively the reflected 
and/or the scattered electromagnetic wave; especially from the sea-surface. 

Finally the shadowing and depolarization effects will be introduced. These effects 
are also significant in determining the reflected and/or scattered field of an electromag- 
netic propagating wave. 

B. SPHERICAL SURFACE AND REFRACTION EFFECTS 

Electromagnetic wave reflectivity becomes more complex when the surface is not 
flat. The reflected wave is decreased in amplitude; also the phase difference between 
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propagation over a convex surface. The difference between the plane and spherical 
surface is that in spherical surface the height of the transmit and receive antennas is 
increased in order to achieve equivalent line of sight. Figure 3.2 illustrates the reduction 
in the reflected wave’s amplitude towards a specific direction due to the divergence of 
the waves into a larger area, compared with a flat surface. [Ref. 2] 

Realizing that the propagation of the electromagnetic wave is in the atmosphere, 
one must consider the refraction effect. This effect causes the waves to bend during 
propagation rather than propagating in straight lines, usually resulting in an extension 
of the maximum range of propagation. It is convenient to analyze the effect of the 
phenomena together with the earth’s curvature. Detailed discussion on atmospheric 
refraction is carried out in Chapter V. 

1. Horizon Range 

The reduced heights {h\,h' 2 ) of the radar and target for a plain surface are 
calculated by: 



h\ = h 1 - 



S- 

2kr„' 



(3.1) 



and 

K = h 2 - (3.2) 

2kr e 

where d\ is the curved radar horizon when h\ is zero, d 2 is the curved distance from 
the reflecting point to the receiving antenna and hi and h 2 are the actual heights of 
the transmitting and receiving antennas respectively (see Figure 3.1). The average 
effect of standard atmospheric refraction is interpreted by increasing the earth radius 
r e by the factor k = 4/3. (For more details see Chapter V.) The differences in heights 
of the radar antenna and the target, due to refraction and earth curvature, are given 
by [Ref. 2]: 
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Figure 3.1: Spherical Surface Geometry [Ref. 2]. 
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Figure 3.2: Divergence by a Convex Surface [Ref. 2]. 
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(3.3) 



d\ 2 

A h 12 = -¥■ 

i, i 2 



The radar horizon d h is given by: 

d h = + hj + \/2fcr e /i 2 + h 2 . 



(3.4) 



The distance d\ + d 2 is the horizon range of the radar if h 2 is zero. Using 
this fact in Equation 3.3, the horizon range then can be approximated to: 



\j2JcrJi, 



(3.5) 



where h is the antenna height of the radar and kr e is the effective earth radius. For 
standard atmospheric conditions ( k = 4/3) the distance dh has the form of: 

d h = 1.23^/h (ft) in nautical miles. (3.6) 

Furthermore, the relationship between the optical, radar and geometric (prop- 
agating along perfect straight lines) horizons for standard atmosphere is: [Ref. 8] 



dh = 1.07d o = 1.15 d g , 



(3.7) 



where dh,d 0 and d g are the radar, optical, and geometrical horizon ranges. Equation 3.7 
stresses that under standard conditions the radar horizon range is 7% larger than the 
optical horizon and 15% larger than the geometric horizon. In abnormal conditions such 
as ducting, the electromagnetic waves can propagate to longer ranges due to diffraction 
from the earth, combined with refraction from atmospheric layers. [Ref. 2] 



C. COMBINED REFLECTION COEFFICIENT 

As mentioned above the physical principle which allows radar coverage to be 
extended beyond the geometrical horizon is divergence or refraction. The electromag- 
netic wave is diffracted around the earth’s curvature, the range of propagation depends 
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on the frequency. The larger the wavelength (lower frequency) the more the wave is 
propagated, due to refraction. 

This phenomena is applied in VLF (Very Low Frequency) for world-wide com- 
munications. Since radax frequencies have short wavelengths compared with earth 
dimensions, only small amounts of energy will be diffracted, resulting in insufficient 
propagation beyond the radax horizon by divergence. [Ref. 5] 

The combined reflection coefficient of eaxth is determined by three factors: 



1) D , the divergence factor that describes the reduction in reflection due to eaxth ’s 
curvature (see Figure 3.2), 

2) p, the reflection coefficient of a plane smooth surface, and, 

3) R , a factor which defines the surface roughness. 

The combined reflection coefficient is the product of the three factors mentioned 
above [Ref. 2]: 

r 

combined 

reflection ’ = pDR. (3-8) 

coefficient 

The factor p may be calculated from Fresnel’s equations (see next section). The 
factor R is equal to unity if the surface is smooth. If the surface is extremely rough R 
is equal to zero. The factor D varies between 0 and 1, and D approaches unity when: 

1) hi or h 2 (Figure 3.1) approaches zero or k becomes very large (especially over 
water), and, 

2) For radax geometry such that the eaxth cam be assumed flat [Ref. 2]. 

The divergence coefficient D is given by: 



D — 



1 + 



2 d\ d 2 



1 - 1/2 



kr e dh sin 6 



(3.9) 



where 6 is the grazing angle. 
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D. REFLECTION COEFFICIENT OF PLANE SMOOTH EARTH 



Prior to discussing the characteristics of reflection of a plane wave from rough 
surface, one should first consider the reflection characteristics from a smooth surface. 
The relative complex dielectric constant t TC is given by: 

e rc = — — 60j‘A<7, (3.10) 

Co 

where eo is the dielectric constant of free space, A is the wavelength in meters, a is the 
conductivity in mho — m -1 ) and e is the dielectric constant. 

The relative magnetic permeability is given by [Ref. 1]. 



Pr=~, (3.11) 

PO 

where [i is the magnetic permeability of earth and p 0 is that of free space. As mentioned 
in Section C, the reflection coefficient of a plane smooth surface can be calculated by 
Fresnel’s equations [Ref. 2,3]: 



_ sing' - (e rc - cos 2 a ') 1 / 2 _ 

• , / o \l/2 P Ht ? 

sin a + (€ rc — cos 1 ayi* 



(3.12) 



e TC sin a 1 — (e rc — cos 2 a') 1 ! 2 .. 

f V ' t \ ( 2 /M /2 i 

e TC sin a' + (e rc — cos^ a ) l/z 



(3.13) 



where a 1 denotes the incident and reflected angle from a tilted surface (see Figure 3.3), 
and 4>h and 4> T are the horizontal and vertical phase of the reflected wave. The mag- 
nitude and phase of the reflection coefficient depends on the polarization, frequency 
and angle of incidence of the electromagnetic wave; and on the dielectric constant and 
conductivity of the surface. 

The Fresnel Equations 3.12 and 3.13 correspond to two polarizations of the electric 
field. The subscripts H and V stand for horizontal and vertical polarizations respec- 
tively. Fresnel equations are given under the assumption that the wave is propagating 
in free space, and that the dielectric surface is non-magnetic. 
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When a is equal to | (normal incidence), the horizontal and vertical reflection 
coefficients are identical in magnitude but opposite in sign. In this case Equations 3.12 
and 3.13 become: 



r 1 - w 

r " I74P 


(3.14) 


e - e 1/2 

r> c rc c rc r> 

i V = = -iff. 

Crc 1 Crc 


(3.15) 



E. ILLUSTRATIONS OF REFLECTION COEFFICIENT OF A 
SMOOTH SURFACE 

Figures 3.4 to 3.7 illustrate the dependence of the magnitude, p, and the phase, 
(f>; on the grazing angle 6 , wavelength A, and the polarization. In the case of vertical 
polarization there is a sharp decrease in the reflection coefficient due to the increase 
of transmission into the surface when the angle of incidence is close to Brewster angle, 
(ag). At this angle py reaches its minimum. The phase <f>y is [Ref. 2,3]: 

(f> v = — — as- (3.16) 

In Figure 3.6 it is shown that <j>y depends also on the wavelength of electromagnetic 
wave. 

When the electric field is horizontally polarized, only slight variations occur in the 
magnitude and phase of the reflection coefficient (pn and </>// respectively), with respect 
to the incidence angle (See Figures 3.4 and 3.5). The phase (j>u is almost 180° for long 
wavelengths and increases slightly with incidence angle. In many practical applications 
pH and 4>n are approximately unity and 180 degrees ( 7 r radians) respectively, i.e., 
r = -1. [Ref. 1,2] 

For very small incident angles py and <f)y are also close to unity and 7r, respec- 
tively. Such a case occurs for anti-ship sea-skimming missiles guided by a radar and for 
microwave searching radars operating close to the ocean surface. For incidence angles 
larger than 10°, <j>y decays to zero (see Figure 3.5). The characteristics mentioned 
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Figure 3.3: Incident and Reflected Wave by a Smooth Dielectric Surface 
[From Ref. 2]. 
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above for small incidence or grazing angles (< 1°) will play a major role also when the 
surface of incidence is rough. A complex effect of shadowing arises for small incidence 
angles and rough surfaces. This effect will be discussed in the next section. 




Figure 3.4: Reflection Coefficient of Very Smooth Sea Water [Ref. 2], 



F. THE EFFECT OF SURFACE ROUGHNESS ON THE REFLECTION 
OF ELECTROMAGNETIC WAVES 

In Chapter II we discussed the Rayleigh criterion and Fresnel zones for smooth 
surface approximation. In this section we will try to consider some more aspects of 
distinguishing between a rough and a smooth surface. Using Equation 2.4 for Rayleigh’s 
criterion of smooth surface we get [Ref. 2,3]: 



AhsinB 1 
A < 8' 



(3.17) 



From this simplified criterion one would expect the reflection coefficient to be a 
function of wavelength, surface roughness, and grazing angle. This criterion gives only 
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Figure 3.5: Phase of Reflection Coefficient <f> of Very Smooth Sea [Ref. 2]. 




Figure 3.6: Reflection Coefficient for Low Depression Angles of Very Smooth 
Sea [Ref. 2]. 
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Figure 3.7: Phase of Reflection Coefficient for Low Depression Angle of 
Very Smooth Sea [Ref. 2]. 



an approximate indication, deeding only with orders of magnitude. From empirical data 
it had been indicated that for large values of the parameter in Equation 3.17, where: 



Ah sin 6 
A 




(3.18) 



the reflection coefficient is independent of the parameter A - ^*' n - e for a specific surface. 
However, small values of A ^ ing , the reflection coefficient increases with an increase 
of the parameter. Actually there is no sharp distinction between rough and smooth 
surface. [Ref. 4] 

One fundamental result from analyzing the scattering field from a rough surface is 
that it consists of two components, a coherent or specular component and an incoherent 
or diffused component. [Ref. 1] 

The specular reQection is defined as a single reflection of electromagnetic wave 
from a rough surface in the direction of the receiver. Its reflection obeys the optical 
laws end its phase is coherent. In other words, the reflected waves from a rough surface 
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can be considered as specular if they arrive at the receiver with equal phases and have 
relatively small amplitude fluctuations. These reflections are from the Fresnel zones as 
described in Chapter II. 

The diffused scattering is scattered randomly from a much larger area than the 
first Fresnel zone. Its phase is uniformly distributed, but the large fluctuations of its 
amplitude are Rayleigh distributed. [Ref. 1] 

These two components were observed in several experiments. The specular re- 
flected field, when combined with the direct field, will produce at a test receiver a 
periodically varying field. When a diffuse field is added to the direct field, the result 
will be a randomly varying field. [Ref. 1,3] 

Specular reflection and diffuse scattering often occur simultaneously. The diffused 
scattered field could be the reflection from large numbers of uncorrelated scatterers. 
It would be expected that the reflection is totally diffused and independent of phase 
difference A<^. This is not always true, since in the case of perfectly smooth surfaces 
the diffused field is zero [Ref. 4]. Also the specular field many be produced by many 
scatterers since it can be produced by a linear summation of many specular fields. In 
reality, the total field is the sum of both the specular and diffused fields. 

Many attempts were made to model a reliable specular reflection coefficient, 
caused by a Gaussian distributed surface as a function of mean wave height of the 
sea, grazing or incident angle 6 , and the transmitting wavelength A. The improved re- 
flection coefficient p s does not take into account the shadowing effect. The magnitude 
of the specular reflection coefficient is calculated by [Ref. 2]: 



p, = p 0 -DR,, (3.19) 

where R s is the specular scattering coefficient for rough surface, po is the reflection 
coefficient of a plane and perfectly smooth surface, and D is the divergence coefficient 
due to earth’s curvature. Modifications in the reflection coefficient caused by a Gaussian 
distributed rough surface was calculated by Beckman [Ref. 1]. The mean square 
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variations in R, is given by: 



<R]>=e~^\ 



(3.20) 



The phase difference A <f> is calculated by: 

A<f> = 2nAhsinO/\, (3-21) 

where Ah represents the standard deviation of the Gaussian distribution of wave heights 
[Ref. 2]. Equation 3.21 is an analogy of phase difference for Rayleigh criterion (Equa- 
tions 2.2, 2.4). The improvement in Equation 3.21 is that Ah is not just a magnitude 
but defined more correctly. 

Most of the properties of reflection from a smooth surface can be applied to 
specular reflection from a rough surface. The basic difference is that specular reflection 
[Ref. 1]: 

1) fluctuates, and, 

2) is smaller in amplitude. 

These facts are taken into account through the specular scattering coefficient R„. 

Different models for < R 2 > were suggested. In all of these < R 2 > is a function 
of the phase difference A <f>. The additional models of < R% > are [Ref. 1]: 

1) a uniformly distributed surface: 

< R] >= sinc 2 A<£, (3.22) 

2) first Fresnel zone of a uniformly distributed surface: 

< R] >= (1 - A<f>/n)smc 2 A<f>. (3.23) 

By checking these three models (Equations 3.20, 3.22, 3.23) against experimental 
results, the improved Gaussian surface distribution, Equation 3.20, gave the closest 
results to the experimental ones [Ref. 2]. Figure 3.8 illustrates the dependence of 
the mean square of the reflection coefficient, on the phase variance, for the improved 
Gaussian model. 
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From experiments it was seen that there were some variations of < R% > from the 
curve when the abscissa is > 120 mradians. Those variations were caused, probably, by 
diffused scattering; or by the combination of diffused and specular, scattering/reflection 
[Ref. 2]. 

Until now we have discussed, quantitatively, only the specular reflection from 
a rough surface. Now we want to define and analyze the diffused scattering. The 
reflection coefficient of diffuse scattering pd can be written similarly to the specular 
reflection coefficient: 

Pd — PoDRd, (3.24) 

where Rd is the coefficient which defines the irregularities of the surface, and po is the 
reflection coefficient of a smooth surface. The amount of power in the diffused field may 
also be a function of polarization. The values of Rd that were measured in experiments 
is within the range of 0.2-0. 4, while the mean values were between 0.3-0.35. When 
6 = 0 (grazing incidence), A <j) = 0, according to Equation 3.21; and < R? a >= 1. 
Under these conditions, p s = Dp 0 . This means that the specular (coherent) component 
is equal to the value of the field for a smooth surface. As 6 increases, the diffuse 
(incoherent) component increases, until the phase difference approaches the value of 
120 mradians. From this point the diffused component decreases. Also, an increase in 
6 decreases the specular component (see Equations 3.19-3.21). [Ref. 2,4] 

For sea surface, one of the models for forward scattering expresses the electric 
scattered field as a vectorial summation of the specular and diffused fields with zero 
mean. The phase of the specular field is equal to the phase for a smooth sea; its 
magnitude is reduced by the scattering coefficient (Equation 3.19). [Ref 3] 

Measurements were made by Beard [Ref. 12] on the reflected field above sea 
surface, for roughness varying from smooth sea < | till the value of = 

The transmitting wavelengths were 0.86, 4.2, and 5.3 cm. Both polarizations, 
vertical and horizontal, were introduced separately. The measured results indicated 
that [Ref. 2]: 
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Figure 3.8: Dependence of Mean Square of Scattering Coefficient on Phase 
Variance, for Specular Reflection from Rough Surface [Ref. 2]. 
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1) The maximum values of power contained in the diffused field is approximately 

30% of the specular reflection. 

2) For both polarizations, the results did not change the polarization independence. 

Probably the most common distribution of forward scattered field from a sea 

surface is the Ricean distribution; which is a result of a model where the amplitude 
fluctuation of the forward scattered diffused field is Rayleigh distributed, and its phase 
is uniformly distributed. This model has proven to be a good approximation, also, by 
experimental results. It should be noted that for veiy low grazing angles, the ratio of 
the diffused to specular powers is 2-3 time larger than for higher incidence or grazing 
angle. [Ref. 3,6] 

The root mean square (RMS) of the sea wave height Ah cam be deduced by 
observation or from sea state tables. These values are given in Table 3.1, and will be 
used in Chapter IV to calculate the coverage pattern of the radar due to multipath. 

TABLE 3.1: PARAMETERS OF SEA STATE. [REF. 1,7]. 



Sea 

State 


Type of 
Roughness 


Wave 
Height [m] 


RMS Height 
[cm] 


0 


calm 


0 


0 


1 


smooth 


0-0.3 


0-6.5 


2 


slight 


0.3-1 


6.5-21 


3 


moderate 


1-1.5 


21-32 


4 


rough 


1. 5-2.5 


32-54 


5 


very rough 


2.4-4 


54-86 


6 


high 


4-6 


0.86-1.3 



G. SHADOWING EFFECT 

When radar operates at very low grazing angle (< 1°), the scattering from the 
surface will be “partial”. This means that large parts of the electromagnetic wave lie in 
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the geometric shadow of the sea wave crest [Ref. 3] . The trough of the wave is not seen 
directly, though some of the shadowed areas are seen due to refraction illumination. 

However, for grazing angles less than 1° the shadow effect greatly affects the 
backscatter echo. 

Two analytical approaches were developed for the shadowing problem. A statisti- 
cal exercise in geometrical optics and a diffraction problem in which geometrical optics 
shadowing theory is considered. The procedure of the statistical approach involves 
many statistical expressions, but eventually results in a shadowing function which is 
the percentage of the illuminated surface area at a specific grazing angle. At small 
angles it is a linear function of the angle. This approach has not yet been proven as 
representing reliable practical shadowing effects at low grazing angles. [Ref. 6] 

In the geometrical shadowing theory mentioned above, it was found that at ex- 
tremely low grazing angles (< 1°) the sea surface resembles a dark plane with illumi- 
nated “spots”. The size and location of these spots, within the reflecting plane, are a 
function of the wind speed and the grazing angle. This result is closer to the practical 
radar picture, involving such phenomenon as sea spikes, rather than the shadowing 
function picture which presents diffuse, vanishing, scatter. However, this geometrical 
theory should be checked carefully before applying it to practical cases. [Ref. 2,6] 

H. DEPOLARIZATION 

Depolarization is defined as a change in the polarization of an electromagnetic 
wave due to propagation (in a medium other than free space) and including surface 
diffraction, reflection, and scattering. 

In describing quantitatively the depolarization scattered from an object, a so- 
called scattering matrix is derived. The scattering matrix describes completely the 
phase and amplitude of the reflected electromagnetic wave at the receiver. For better 
understanding, suppose that an electromagnetic plane wave is transmitted with its 
electric field in the y direction (see Fig 3.9), and represented by E y . Let the reflected 
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field from an object (or a surface) at the receiver, be represented by E x , and be given 
by [Ref. 8]: 



E x — P VX E V , 



(3.25) 



where p yx represents the reflection coefficient of the incident surface and/or object. 
(The first subscript stands for the transmitted field and the second subscript stands for 
the receiving component.) The reflection coefficient p yx is complex, and includes the 
change in phase and amplitude, of the incident field. 

Since any polarization of a plane wave can be described by two orthogonal com- 
ponents, polarized linearly, a scattering process can be described completely by the 
polarization scattering matrix [after Ref. 8]: 



H -l 

- 




Pxx 


Pyx 




' Ei' 


E r 

v J 




Pry 


Pyy . 




. E l. 



(3.26) 



Figure 3.9 illustrates two linear and orthogonal electric fields, which produce by sum- 
mation a right circular electric field. For a perfect conducting plane surface, p xx should 
be equal to Also, the scattered energy (square electric field) from a surface is 

mainly a function of the transmitted wavelength and the reflecting surface orientation 
and roughness. [Ref. 8] 

In 1963 it was shown by Beckman and Spizzichino [Ref. 1] that, if a plane wave is 
purely polarized (horizontally or vertically), no depolarization will occur in the reflected 
wave in the plane of incidence. If the same incident wave will be scattered out of the 
incident plane , it will be strongly depolarized. These scattered components are a 
result of the irregularities of the incident surface, and related to the diffused scattering. 
[Ref. 1] 
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Figure 3.9: Right Circular Polarization Produced by Horizontal and Verti- 
cal Polarizations [Ref. 8]. 
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According to experiments above sea surface, with near purely polarized wave, the 
received depolarized specular component is very small (at least 25 dB below the incident 
polarization). On the other hand, the diffused scattered field is strongly depolarized. 

From these results one can conclude that for forward scattering from sea surface, 
depolarization, for linear polarization is small. This fact can also be verified when 
looking at echoes reflected from targets, especially on a smooth sea. These echoes 
might fluctuate sharply due to interference of the direct and the reflected wave. The 
interference can occur only if the reflected field has largely the same polarization and 
magnitude as to the direct field. It is shown in Ref. 1 that very little depolarization 
takes place in the direction of specular rays. Therefore the interference effect is addi- 
tional evidence that depolarization, due to reflection from sea surface is small [Ref. 2]. 
This phenomena will be discussed and analyzed in Chapter IV. 
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IV. LOBING PHENOMENA ANALYSIS 



A. INTRODUCTION 

In Chapters II and III we have described the characteristics of earth’s surface and 
a little about air refraction (for more detail see Chapter V), which dramatically affects 
electromagnetic propagation. 

In this chapter we will discuss and analyze the influence of these characteristics 
on a Naval search radar (shipborne, shore installation and low-flying aircraft). 

We can assume that the detection range of a very low altitude target will be 
limited due to multipath propagation and earth’s curvature [Ref. 11]. 

In this chapter we will discuss and analyze the lobing phenomena, caused by 
multipath, based on the parameters described in Chapter II and III. We will consider 
specifically the sea surface, and the detection conditions of a Naval search radar, which 
usually views the sea at extremely low grazing angles (< 1°). The covering pattern 
of detection for these radars is strongly affected by the multipath of the propagating 
electromagnetic wave. The multipath causes constructive and destructive interference 
due to earth reflection. This effect, together with media and radar parameters, will be 
computed to illustrate predicted detection patterns of the radar. 

We will complete the analysis with some practical conclusions and recommenda- 
tions for some control on the pattern coverage for a Naval radar. 

B. MULTIPATH FADING PHENOMENA 

The main sources of multipath fading are reflections from the earth’s surface and 
from refracting layers in the atmosphere. In the commonly known multipath, there are 
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two such paths-a direct path of the transmitting wave from the radar to the target and 
an indirect path from the radar to the target via reflection. 

The geometry of this scenario is shown in Figure 3.1. From this figure it is seen 
that the direct and indirect wave have traveled different paths. If these waves have the 
same order of magnitude they can strongly interfere. They will have phase differences 
caused by path differences. A quantitative analysis of the lobing phenomena was done 
by Blake and is described herein [after Ref. 9]. 

Destructive interference will occur when the absolute value of the phase difference 
is between | + 27rn and + 27 rn radians (n is integer). When the phase difference is 
it radians (180°) and the amplitudes of the direct and indirect field are almost the same, 
a complete cancellation will occur at the receiving point. Constructive interference will 
occur if the absolute value of the phase difference is between 0 + 27rn and 7r/2 + 27rn. 
When there is no phase difference (<f> = 0), the constructive interference will be of 
maximum value due to arithmetic summation of the amplitudes. 

This phenomena can be seen, if we look at the field at a fixed height, as a function 
of the range. We will cross through maxima and minima of the field pattern of the 
radar, caused by interference. 

The goal of this discussion will be to derive procedures and equations in order to 
predict radar range performance under multipath conditions. We will derive the lobing 
factor F ; considering the flatness of the earth, sea roughness and polarization. This 
was discussed in Chapters II and III. For simplification, we will assume that the target 
is far enough so that the usual fax-field approximations are valid. 
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C. THE LOBING FACTOR 



The analytical treatment of the interference phenomena has to be in the complex 
plane, since we are dealing with amplitudes and phases of the waves. Basically, the 
phase difference of the direct and reflected wave is dominated by three factors: 



1) the geometric path length difference 8 (see Equations 2.6, 2.7), 

2) the phase difference of the reflected wave due to the reflection process, and, 

3) phase difference due to radiation of the antenna in direct and indirect directions. 
The amplitude difference can be caused by the following reasons: 

1) The Naval radar with narrow beam in the horizontal plane is usually not uniform 
in the vertical plane which can cause different field intensities in the direct and 
indirect waves. 



2) Both waves propagate spherically, resulting in an inverse square law power reduc- 
tion. Since the path of the waves is different they will be attenuated or spread 
differently. 

3) The reflected wave is usually attenuated because the sea surface reflection coef- 
ficient is less than unity. 

1. Applications 

The spreading loss is usually neglected for both waves due to the small path 
difference 8 between the direct and reflected wave. The major effect of 8 is on the phase 
difference. Most of the Naval surface surveillance radars Eire pulsed radars. The effect 
of interference can occur only if 8 < ct, where r is the pulse width and c is the speed 
of light in free space. For surface and low flying targets and moderate height of radar 
antennas, 8 <C cr [Ref. 5,9]. This will result in susceptibility to the multipath effect. 
For radars which use pulse compression techniques the criterion of 8 < ct is for the 
compressed pulse. 

2. Combined Phase Difference o, and Detecting Range R 

The magnitudes of the reflection coefficient p, and the pattern factor of the 
antennas /(7), axe major effects on the amplitude of the reflecting wave. The elevation 
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angle between the radar antenna and the target/receiver is 7 (see Figure 4.1). The 
ratio of the electric intensity in the direction of the reflected ray to the intensity of 
the maximum beam direction (boresight) is /( j r ), while f(id) is the same ratio in the 
direction of the direct ray. Both 1 /( 7 ) | and p are bounded in the range of 0 and 1 . 

Now we can derive a phasor expression for the direct and indirect wave fields. 
The direct wave field can be written as: 



E, = fMEo, (4.1) 

where Eo is the phasor value of the free space field in the beam maximum at the target 
range Rd. Similarly, the reflected wave field is: 

E r = f(y r )Eo-pe-> a , (4.2) 

where a is the total phase difference caused by three factors mentioned in Section C of 
this chapter, and p is the magnitude of the reflection coefficient. The combined field at 
the target is: 

E = E r + E d , (4.3) 

Substituting Equations 4.1 and 4.2 into 4.3 yields: 

E = E 0 {fh d ) + (4.4) 

the lobing factor F is defined as: 

F = \E/E 0 \. (4.5) 



Dividing Equation 4.5 by E 0 yields: 



F = \fM + pfMe-n 



(4.6) 



If we assume that o includes the phase difference (3 r — /?<* of the pattern factor, 
then by factoring out f d and using Euler’s theorem, the lobing factor F is written as: 



F = /j|\/l + x 2 + 2xcoso|, 



(4.7) 



38 



where 



(4.8) 



X = 



Pfr 

fd' 



f T and fd are scalar magnitudes of the direct and reflected pattern factors. 

We can represent the pattern factor f(yd) and f(lr) in the form of: 



f(-u) = (4.9) 

fM = (4.10) 



where /3d and P T are phase angles and fd and f T are the magnitudes of the pattern 
factors. 

As mentioned, ot is the total phase difference between the direct and reflected 
waves and can be expressed as: 

ttS 

a = — + <f> + fir - Pd, (4.11) 



where the first term in a is caused by path difference, the second term (<^) is caused 
by the reflection process, and the third term ( P T — Pd) is the phase difference of the 
pattern factors. 

The detection range in presence of multipath and atmospheric loss L can be 
calculated by: 



R = 



RoF 
1 'A/*' 



(4.12) 



where Rq is the free space detection range. For analyzing the lobing effect we will 
propose that no atmospheric loss occurs, and Equation 4.8 is simplified to: 



R = FRq. (4.13) 

Equation 4.13 is a realistic approximation for radar frequencies, since the 
lobing factor is dominant when multipath occurs and affects the detection range more 
than the atmospheric losses. Since 0 < F < 2, it is seen that the range Rd under 
multipath conditions is sometimes greater than the free space range, and sometimes 
smaller than the free space range. 
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D. PATH DIFFERENCE DERIVATION 

The actual antenna pattern located at A in Figure 4.1 can be found by treating 
it as a two source interferometer caused by the real antenna and its image from the 
ground at C [Ref. 5]. The range Rj is the flat earth approximation range. If we assume 
that Rd h[ then the difference between the direct path AE and the indirect path 
ADE is: 

8 = 2h[ sin 8, (4.14) 

Recall from Chapter III, Section A, h\ is the height of the radar antenna and 8 is the 
incident (grazing) and reflected angle of the indirect ray. 

For small grazing angles, which are typical of sea surface search radars, sin 8 can 
be expressed as: 



sinf? = 



h\ + h' 2 

Rd ’ 



(4.15) 



where h[ is the height of the target. By assuming that h 2 h\ and knowing that 
Rd h[ the further approximation can be made: 



a K 
sin 8 ~ — . 

Rd 



(4.16) 



By substituting Equation 4.16 into 4.14, the path difference can be approximated by: 



2 h[ ■ h' 2 
Rd 



(4.17) 



By multiplying Equation 4.17 with the wave number, k = the result is the approx- 
imate but practical phase difference A <^>, for a plane reflecting surface: 

27r 2 h\ • h 2 



A<f> 



A 



Rd 



(4.18) 



The elevation angle jd, and the depression angle y r , can be calculated by geometric 
considerations: 

h' — h' 

(4.19) 

•fid 

and 



7 r = sin 



R d ’ 

-i ^2 *b h\ 
R\ + R 2 



(4.20) 
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Figure 4.1: Multipath Geometry over a Plane Reflecting Surface [Ref. 9]. 
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Again by using the assumptions h' 2 h\ and Rj h\ we conclude that : 



Id = 7r = 0 . 



(4.21) 



By using the identity of Equation 4.21 in Equation 4.14 there results: 

6 = 2 h[ sin(7d). 



(4.22) 



This geometric solution of the multipath effect is not complete since we have to 
consider also the reflection coefficient which is a function of the grazing angle 9. 

E. PRACTICAL LOBE FACTOR DERIVATION 

The solution for the lobing factor is achieved by substituting Equation 4.17 into 
Equation 4.11; this result for a is substituted into Equation 4.7 together with the 
magnitudes of the pattern factors and the reflection coefficient. [Ref. 3,9] 

A special, practical, and simplified result is obtained for a perfectly smooth sea 
and horizontal polarization. For this case, Figures 3.6 and 3.7 indicate that p ~ 1, and 
<f> ~ 7r rad, independent of grazing angle. 

For surface search and navigation radars the main beam is narrow and directed 
at zero elevation angle, and the vertical beam is broad. Under these conditions (3^ = (3 T 
and fd = /r = 1 and 

9 - 7 rf\ 

(4.23) 



2tt6 

a = — f- 7 r. 



Substituting this simplification into Equation 4.7 yields: 



F = |^2 + 2 cos(27r6/'A + 7r ) |. 
By trigonometric manipulation: 

F = 2| sin(7r£/A)|. 



(4.24) 



(4.25) 



By substituting Equation 4.21 into 4.24 we obtain, finally, an especially useful result 
for the lobing factor, F : 



F = 2 



sin |(27r/ii sin 7 d)/A] |. 



(4.26) 
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F. DISCUSSION 



Equation 4.26 allows us to predict the relative lobing pattern without knowing 
the range and altitude of the target. The lobing factor F depends only on the elevation 
angle 7 <*. Figure 4.2 is a plot of the lobing factor F showing the vertical plane coverage 
affected by typical multipath from a smooth sea. The antenna height is 10 meters and 
the transmitting frequency is 100 MHz. 

Note that the maximum range Rd is twice the free space range Rq due to con- 
structive interference since the maximum value of F is 2. On the other hand, when the 
target is located at other elevation angles than the maximum beam value, the detec- 
tion range can be less than the maximum free space range. When the target is in the 
direction of nulls, no detection will occur. 

As seen in Figure 4.1, variations of detection range relative to free space are 
bounded between the two extremes, 0 and 2, and are a function of the elevation angle. 
At a lobe maximum, we have [Ref. 5,9]: 

7max(n) = sin -1 A n = 0,l,2... . (4.27) 

The minimum occurs when the same argument is equal to multiples of 27rn. The 
elevation angles at which minima occurs is given by: 

nA 

7min(n) = sin -1 — n = 0,1,2... . (4.28) 

Equations 4.27 and 4.28 indicate that the maxima and minima angles depends 
only on the wavelength and antenna height. This is true only under the following 
assumptions [Ref. 9]: 

1) The pattern factor phase difference /?<* — /? r is constant at all elevation angles. 

2) The phase of the reflection coefficient is also constant. 

3) The reflecting surface is plane and not curved. 

These assumptions are good approximations for horizontal polarization and an 
antenna pattern with zero elevation beam. If we take into account also the reflection 
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Figure 4.2: Vertical Lobe Pattern Caused by a Plane Reflecting Surface. 
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coefficient and antenna pattern factors, it affects the envelope of the covering factor, 
as we will see in Section G of this chapter. 

The free space radar equation for the received signal should be modified by the 
lobing factor. Due to free space loss the received power P r of a radar (two-way loss) is 
proportional to: 

P, oc k. (4.29) 

Doubling the free space detection range Ro caused by constructive interference is equiv- 
alent to increase the transmitting power or the sensitivity of a radar 16 times or equiv- 
alently by 12 dB. 

Equation 4.28 yields that the first minima occurs at zero elevation angle. This is 
not true for a curved surface in which the first minimum is not zero and its elevation 
angle is not zero. Near the horizon when the grazing angle approaches zero, the lobing 
factor is very low causing a very low field strength at this region. The returned signal 
is dominated by the spherical nature of the wavefront, resulting in scattering which 
can be taken into account by the divergence factor D [Ref. 3]. This is the reason why 
targets at very low grazing angles (low altitude), much under the first maximum lobe 
angle are difficult to detect. Typical targets with low grazing angles with respect to 
shipborne radars are anti-ship sea skimming missiles. This is the reason for the effort 
being made trying to predict and control the location and shape of the first interfering 
lobe. 

The flat earth equations that were presented in this chapter are good approxima- 
tions for relatively low antenna height such as shipborne searching radars where the 
reflecting point is close to the radar. 

G. LOBING ANALYSIS FOR TYPICAL CASES 

We want to now analyze plots of detection patterns affected by multipath for 
radars operating above sea surface. In Figure 4.2 we saw a detection contour of a 
low frequency radar (100 MHz) with horizontal polarization and antenna height of 10 
meters operating above a perfectly smooth surface. Figure 4.3 is the detection contour 
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of the same radar except that the height of the antenna is doubled to 20 meters. This 
has an effect of lowering the angle of the first lobe and doubling the number of lobes 
nulls. 

Figure 4.4 is a detection plot of a radar under the same conditions as Figure 4.3 
except the frequency is now 10 times higher (1000 MHz). This effect also lowers the first 
beam and reduces considerably the uncovered space between the lobes. The number 
of lobes and nulls are increased by the ratio of frequencies (10 times). Also, according 
to conservation of energy, the lobes are 10 times narrower. One can conclude that 
the “penalty” for lowering the first lobe by increasing the antenna height and/or the 
frequency is that the number of nulls is increased proportionally (no received echo in null 
angles) and that the lobes are narrowed proportionally. For shipborne and shorebased 
surface search radars, increasing the frequency and the antenna height is recommended 
for improved detection of targets with low grazing angles like small surface combatants 
and low flying aircraft/missiles. 

If the polarization of the transmitting wave is vertical, the covering pattern will 
be reduced. This is particularly true for the lower lobes, due to reduced values of the 
reflection coefficient p, at low angles (see Figure 3.4). Also, the phase difference is not 
180° as it was in the case of horizontal polarization (see Figure 3.5), resulting in higher 
angles of maxima and minima. [Ref. 9] 

In Figure 4.5 we include also the effect of the antenna pattern of a rectangular 
antenna with its electric field uniformly distributed. The pattern is derived from the 
Fourier transform of the rectangular uniform distribution resulting in a shape. By 
comparing with Figure 4.3, we see that this has the effect of modulating the basic 
pattern. Also the higher order lobes have a reduced coverage. 

We now investigate the influence of the scattering coefficient Rs, considered in 
Chapter III, on the radar range. The effect of this parameter is to reduce the detection 
range progressively at higher grazing angles. This is because < R] > decreases expo- 
nentially with increasing grazing/elevation angle 9. Recall from Equations 3.20 and 
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Figure 4.3: Vertical Lobe Pattern for the Same Parameters as Figure 4.2 
Except Antenna Height Increased to 20 m. 
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Figure 4.4: Vertical Lobe Pattern for the Same Parameters as Figure 4.2 
Except the Frequency is Ten Times Higher (1000 MHz). 
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3.21 that the mean square variation in Rs is given by: 

< R] >= (4.30) 

Figure 4.7 is a plot with same parameters as in Figure 4.6 but for a higher 
roughness of sea surface. By increasing the roughness of the sea-mean wave height 
to 0.6 meters instead of 0.25 meters, the entire pattern is considerably reduced due to 
an increase of the diffuse scattering over specular coefficient. 

In Figure 4.8 we want to show the detection contour of a radar with appreciably 
higher frequency (1000 MHz instead of 100 MHz). The elevation angle is also reduced 
since the beamwidth is proportional to where D is the physical dimension of the 
antenna. The plot is carried out for a perfectly smooth sea and horizontal polarization. 
We can see the effect of the symmetrical antenna pattern on the detection contour. 
The effects of higher frequency are [Ref. 5,9]: 

1) the lobes are finer, 

2) more closely spaced, and, 

3) the first lobe is at lower elevation angle. 

The effect of sea roughness at high frequencies is slightly different from that at 
low frequencies. Figure 4.9 uses the same radar as in Figure 4.8 except mean sea wave 
height is changed from 0 to 0.6 meters. At high frequency the roughness of sea has 
a stronger effect especially on the higher angle lobes. An additional effect which is 
not included in the program of these plots is that nulls are partially filled in this case, 
especially in the side lobes. The multipath effect disappears due to a considerable 
decrease in the specular reflection coefficient p t at high angles [Ref. 10]. 

If the signal is spread over a band of frequencies, the nulls may occur around 
one frequency and not in another frequency, resulting in a decrease in fading time. By 
“spreading” the multipath effect over a wide band (whitening it, in a sense) we can 
reduce the effect of no detection at nulls which occur for radars operating at a single 
frequency. 
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Figure 4.5: Detection Pattern of Rectangular Antenna with Uniform Field 



Distribution. 
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Figure 4.6: Detection Pattern Same Parameters as Figure 4.5 except In- 
cluding the Scattering Coefficient Effect. 
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Figure 4.7: Detection Pattern Showing the Effect of Sea Roughness. 
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Also at 1000 MHz if we change the polarizations to vertical, the lobes magnitudes 
would reduce at elevation angles higher than 2° as can be inferred from the pv curve 
in Figure 3.4. 

Some typical Naval search radars operate in the S band (2-4 GHz). Figure 4.10 
illustrates this practical detection pattern at ideal conditions for lobing phenomena 
over a perfectly smooth sea. 

A very useful way to analyze the multipath effects on the detection capability of 
a radar is a plot of the received signal versus range. This method is especially used to 
analyze targets flying at constant altitude towards a radar. This situation can simulate 
the detectability of an anti-ship missile approaching a ship. 

This plot is more convenient for analysis of detection capability of targets in the 
first lobe (very low altitude targets) within the horizon range. In this plot —93 dBm 
is the minimum detectable signal (MDS). Signals below this level will not be detected. 
We see that the “dips” in echo signal due to multipath are narrower and more frequent 
as range decreases. This effect has been experimentally verified. 

We will only mention here that surface reflection causes errors in height measure- 
ments. This is because the reflected wave can be considered by the radar as an echo 
signal from a target located at the same range as the actual target but at the opposite 
height, i.e., the image of the actual target. When constructive interference occurs, the 
radar will measure the mean value of the elevation angles of both the actual and the 
image targets, resulting in a severe error in height measurement. 
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Figure 4.8: Detection Pattern for Higher Frequency (1000 MHz) and 

Narrower Vertical Beam Width. 
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Figure 4.9: Effect of Sea Roughness on a 1000 MHz Detection Pattern of a 
Radar. 
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Figure 4.10: Detection Pattern under Lobing Effect for S Band Naval 

Surface Searching Radar. 
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Figure 4.11: Received Signal Relative to Minimum Detectable Signal, for 
Multipath and Free Space Conditions. 
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H. SUMMARY AND CONSEQUENCES 



We have seen in this chapter that the specularly reflected field may, if it is com- 
parable in amplitude to the direct propagating field, produce deep fades in the total 
field through interference. The diffused field has much less effect on the fading when it 
occurs. 

Atmospheric refractive layers can cause many different path rays. In this chapter 
we considered only the interference caused by the direct ray reflected from sea surface, 
with no effect of the atmospheric refractive layer. 

The main parameters affecting the detection pattern of a radar are: 

1) Antenna height : Increasing the antenna height increases the number of lobes 
with less uncovered spaces in the detection pattern. The first beam is lowered 
resulting in better detection of low altitude targets. 

2) Transmitting frequency : Higher propagating frequency will produce a covering 
pattern with finer and more lobes. The lowest lobe will be considerably lower. 

3) Polarization : The lobing phenomenon is more prominent for the horizontal po- 
larization since the magnitude and phase of the reflected waves are more closely 
correlated to those of the incident wave. For vertical polarization the reflection 
coefficient (magnitude and phase) depends strongly on the grazing angle and 
decreases at higher angles and frequencies, resulting in a decreased detection 
pattern. The effect is stronger at higher frequencies (above 1000 MHz). 

4) How to decrease the lobing phenomena ? For some operational applications, it is 
possible to overcome the multipath effect. For example, in the detection of low 
altitude targets and in active jamming, fading is not desirable. Some means of 
reducing the fading phenomenon are: 

ai)se of spread spectrum techniques (such as frequency agility or hopping), result- 
ing in a spread of fading time over the spectrum. 

b^ansmit vertical polarization at frequencies above 1000 MHz. At very high 
frequencies (35 GHz), the multipath effect is minimal. 

The above parameters can be easily controlled by the radar designer. 
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V. ANOMALOUS ATMOSPHERIC 
PROPAGATION 



A. INTRODUCTION 

In the previous chapters, we considered and discussed mainly the characteristics 
and effects of earth’s surface, i.e., sea surface, assuming normal atmospheric conditions. 

In this chapter we will consider propagation under abnormal conditions called 
ducting or anomalous propagation caused by rapid decrease of the refraction index, 
especially in warm areas such as the Mediterranean Sea. This phenomenon considerably 
extends the detection range of a radar within the ducting layer compared to horizontal 
free space range. On the other hand, from conservation of energy, the increase range in 
one path will cause a decreased range in another path, resulting in holes of detection 
outside the ducting layer. 

We will consider mainly the surface or evaporation duct because it is the most 
likely one above a sea surface. Refraction losses are also discussed. We will end with 
some conclusions for radars operating in a duct. 

B. TROPOSPHERIC REFRACTION 

Electromagnetic waves propagating in the troposphere do not travel in straight 
lines; instead they tend to bend, due to refractivity of the troposphere, which is related 
to its dielectric constant. This dielectric constant is a function of temperature, pressure, 
and water vapor content of the air. [Ref. 9] 

The two main effects of this phenomenon on radar performance are: 

1) An increase in the maximum range of propagation is increased, and, 

2) errors in measurements of elevation angle. 
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Figure 5.1 illustrates a simplified scheme of these effects. 

1. Refraction Effect 

Refraction or bending of electromagnetic waves in the atmosphere is caused 
by the variation of the velocity of propagation with height. 

The index of refraction is the ratio of the velocity of propagation in free 
space to that in an actual medium. For microwave frequencies propagating in the 
troposphere, which contains water vapor, the refractivity can be calculated by [Ref. 5]: 



N = (n - 1)10 6 = 



77.6p 
— + 



3.73 x 10 5 e 

j , 2 



5 



(5.1) 



where p is the barometric pressure in mbars, e is the partial pressure for water vapor 
in mbars and T is the absolute temperature in degrees I\. The refractivity N is used 
instead of the index of refraction n because it is a more convenient unit when dealing 
with propagation. The index of refraction normally decreases with height since the 
barometric pressure p and the water vapor content e decrease much faster with height 
than the temperature T. 

The index of refraction is also given by: 



c 

n = — , 



(5.2) 



where c is the velocity of light in a vacuum (3 X 10 8 m/sec) and v g is the velocity 
of the propagating wave in the actual medium. Near the surface of earth the index 
of refraction is 1.0003. Since n decreases with altitude the velocity of propagation 
increases with altitude. According to Snell’s law the propagating wave tends to bend 
towards the lower velocity. The result is an increase in propagating range. 

Since the refraction index n, varies with height, it can be shown from Snell’s 
law for spherical geometry that the factor k is varying with height and has the form of 
[Ref. 3,5]: 



1 + r e (dn/dh ) ’ 



(5.3) 
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Figure 5.1: (a) Extension of Radar Horizon; (b) Elevation Angular Error 
[Ref. 5]. 
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Figure 5.2: (a) Bent Beam due to Refraction; (b) Straight Line Propagation 
with Effective Earth Radius 4/3 r e [Ref. 5]. 
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where dn/dh is the gradient of the index of refraction with height. 

In Chapter III we considered the refraction effect by increasing the earth’s 
radius by a constant factor of 4/3, to obtain the effective earth radius for straight line 
propagation. This is only an average radius and should not be used for anything other 
than general calculations. It gives a rough idea of refractivity under normal atmospheric 
conditions. Figure 5.2 shows the replacement of eaxth radius by an effective radius of 
4/3 r e assuming a homogeneous atmosphere and constant gradient of index of refraction 
with height. 

The variation of index of refraction with height was found to be an exponen- 
tial function rather than a linear function as assumed for k = 4/3. The exponential 
model for refractivity is given by [after Ref. 5]: 

N, = N e exp [~C e (h T - h R )] , (5.4) 

where N e is the refractivity at the surface of the earth, C e is a constant, h R is the radar 
antenna height and hx is the target altitude. Use of this refractivity model results in 
smaller errors in elevation angle measurements of the target. 

C. REFRACTION LOSS 

At normal refraction conditions, when the vertical gradient of refractivity is con- 
stant = const.), the main effect will be an increased elevation error as measured 
by the radar. However, another effect, which is considered by Weil 1973 [Ref. 9], is a 
divergence loss of signal strength that is caused by the refractions of rays from atmo- 
spheric inhomogenities. This loss is also known as Lens Effect Loss since the effects 
caused by the atmosphere are similar to the lens effect. 

To illustrate this effect consider a source radiating four rays A,B,C,D as in 
Figure 5.3. Without any absorption or refraction disturbances at the distance R, these 
four rays will be spread out through an area equal to SR 2 , where 6 is the angular 
separation in radians between the rays. Now assume that the rays are propagating 
under refractive conditions which will cause them to bend rather than to propagate 
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in straight lines. The rays C and D have a slightly smaller angle with respect to the 
vertical gradient of refraction, thereby causing them to bend slightly more than rays A 
and B. This will result in a larger spreading area of these rays, causing a decrease in 
power density. [Ref. 9] 




Figure 5.3: Ray Pattern for Refractive Loss Analysis [Ref. 9], 

This refractive or lens effect loss, can be calculated by the ratio of wave front 

areas: 

L' = j, (5.5) 

where A! is the wavefront area under refractive conditions and .4 is the free space 
wavefront area. 

Figure 5.4 shows the lens effect loss in dB as a function of the range of propagation 
with the elevation angle as a parameter. The conclusion from this graph is that, 
although the lens effect is not large, it is important in the case of long range detection 
radars trying to detect targets at low altitudes (low angles). At angles higher than 2°, 
the lens effect loss is seen to be negligible. 



G4 




Figure 5.4: Lens Effect Loss as a Function of Elevation Angle vs. Range 
[Ref. 4]. 

D. PROPAGATION UNDER ANOMALOUS CONDITIONS 

Theoretically, under ideal conditions with no attenuation and absorption, the 
vertical gradient of index of refraction could have a high value causing rays to curve 
similar to the earth curvature resulting in continued propagation around the earth. 
The effective radius will be almost infinite if jj- = —1.57 • 10 -7 m -1 allowing horizontal 
rays to propagate around the earth s curvature with no horizon limits. This condition 
is called “ducting,” in which the electromagnetic ray is “trapped.” These conditions 
can occur only within an atmospheric layer with a finite thickness. 

Values of vertical refractive gradient for normal and abnormal propagation are 
known to be [Ref. 5]: 
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1) For normal atmospheric refraction 0 < < —0.787 • 10 7 m 1 or 1 < K < 2. 

2) For superefraction conditions -0.787- 10 -7 < < —1.57 ■ 10 - 7m -1 or 

2 < K < oo. 

Ducting occurs when gjr > -1.57- 10 -7 m -1 . This phenomenon can be compared 
to a leaky waveguide. There are places along the layer in which the index of refraction 
increases with altitude instead of decreasing, causing the ray to bend upwards instead 
of downwards. In spite of the “leakage,” ducting usually increases tremendously the 
horizontal detection range of radars. 

Once duct conditions are fulfilled, the refractive index within the duct is greater 
than the threshold for ducting, resulting in the ray curving down in a curvature larger 
than that of earth’s surface. If the surface is water, some specular reflection will occur. 
Since the ray is trapped, this process will proceed continuously within the duct as 
shown in Figure 5.5. It is shown in the figure that there is a critical elevation angle 
at which the ray propagates within the ducting layer. Beyond this elevation angle, the 
rays will propagate out of the layer. This critical angle is usually less than 0.5° meaning 
that the rays are almost parallel to the surface. [Ref 9] 

This type of duct with a low layer above surface is called a surface duct, or 
evaporation duct, if it occurs above water surface. A duct which occurs at higher 
layers is called an elevated duct. [Ref. 5,9,13] 

The evaporation duct usually occurs within a few tens to almost a hundred feet 
above the sea surface. The duct effect can be interpreted as “forcing” the wave to 
propagate in one dimension (horizontal) rather them in the two dimensions under nor- 
mal conditions (horizontal and vertical). This results in power density spreading as the 
inverse first power of the range (1 /R) instead of l/R 2 . In two way propagation radar, 
the receiving energy will be proportional to l/R 2 under ducting conditions instead of 
the l/R 4 in free space. The result is that for a radar and a target within the ducting 
layer, the detection range will be much beyond the “horizon” range of free space (see 
Fig. 5,6). [Ref. 9, 13] 



66 




t 

LU 

© 

Z 

< 

tr 



Figure 5.5: Ray Propagation Pattern under Ducting Conditions [Ref. 9]. 
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Though there are losses caused by absorption and leakage, the ducting phenomena 
has an effect of increasing the detection range within the layer. This extended detection 
range occurs because more energy is directed into the ducting layer by refraction, while 
less energy will propagate in other directions causing “holes” or poor detection in areas 
outside of the ducting layer. We can conclude that a surface or evaporation duct favors 
detection of surface or sea skimming targets, and does not favor the detection of targets 
with elevation angles greater than 1°. 

E. SURFACE DUCT CHARACTERISTICS AND ATMOSPHERIC RE- 
LATIONSHIP 

As mentioned in Chapter V, Section A.l, the refraction index decreases with 
height in normal atmospheric conditions. Under certain conditions, especially in sub- 
tropical regions, there is a warm layer over a colder layer which means that temper- 
ature increases with height instead of decreasing with height. This, together with 
normal rapid decrease of pressure and evaporation, results in a more rapid decrease 
of refraction index than under normal conditions. Such ducts are common in tropical 
and semitropical areas over cold sea water with a maximum height of at most, a few 
hundred feet. [Ref. 3,9] 

Another type of duct, the evaporation duct, is due to evaporation of moisture 
from the sea surface. This duct is usually lower, 20-60 ft on the average, and is more 
reliable than the previously discussed duct. This type of duct is common in the warm 
areas of the globe such as the Mediterranean Sea, for instance. 

The maximum wavelength of an electromagnetic wave able to propagate within 
a duct depends on its thickness d. A duct could be thought of as a high pass filter 
allowing only the higher frequencies (shorter wavelength) to propagate abnormally. 
The maximum wavelength that can propagate in evaporated duct with thickness d is 
given by [Ref 5]: 

A™=2. 5 (-|)‘V. (5.6) 
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Figure 5.6: Relative Received Signal Under Duct and Free Space Condi- 
tions. 
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For an X band radar (A = 3 cm) the required thickness for ducting to occur, would 
be 10 m. Since ducting layers are not very thick, ducts will usually occur for higher 
microwave frequencies. 

Compared with waveguides, the “cut-off wavelength” for a duct is not so strict. 
Modes with “harmonic wavelength” can also be propagated by ducts. [Ref. 5] 

F. CONCLUSIONS OF ANOMALOUS PROPAGATION 

Finally the following points should be considered when designing and/or operating 
a radar in anomalous propagating conditions: 

1) Evaporation duct, when used with the appropriate height of antenna and wave- 
length, can extend the detection for surface or low flying targets considerably 
beyond the expected range with normal atmospheric conditions. 

2) The extended horizontal detection range will cause reduced range or holes in 
detection at other elevated angles. This follows from conservation of energy. 
This can cause severe coverage problems for air surveillance and early warning 
radars. 

3) Usually radar is designed under the assumption that clutter is likely only in the 
range close to the radar. In ducts, clutter is likely also at longer ranges, resulting 
in a signal-to-clutter problem for radars without techniques to cancel it (such as 
pulse staggering against second time around clutter.) 
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VI. SUMMARY AND CONCLUSIONS 



A. SUMMARY 

In this study we have tried to quantify the propagation characteristics of an 
electromagnetic wave through material and media boundaries. Random irregularities 
of the sea surface, the so-called sea roughness, have been investigated and included in 
the specular and diffused reflection coefficients. The curvature of the earth’s surface 
is taken into account by way of an effective radius. These parameters were used to 
compute the detection performance contours of a typical microwave search radar. The 
multipath effect appears as a lobing factor on the detection range. Further, the effects of 
anomalous propagation through the atmosphere are also studied. The work illustrates 
that the propagating medium has a significant effect on the detection range. 

B. CONCLUSIONS AND RECOMMENDATIONS 

The results of this study show that the reflected field from the sea surface interferes 
with the direct field, and depends on wavelength, polarization, sea state, angles of 
incidence, and the electric properties of the sea surface. 

By appropriate design and operation of the radar, we can have some control on 
the detection pattern. 

Frequency and/or spatial diversity, of the transmitting frequency and antenna 
heights respectively, were proved to be very effective in “filling up” dips in the receiving 
signal caused by multipath. Also, increasing the frequency and antenna height lowers 
the first lobe of the detection contour, for better detection of targets at low grazing 
angles. 
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In areas such as the Mediterranean and the Red Seas, abnormal atmospheric 
conditions for propagation are very common for most parts of the year. In order to 
gain the benefit of the ducting phenomena one must predict the thickness of the ducting 
layer. 

By using appropriate antenna height and wavelength, detection of targets within 
the ducting layer are maximized. Further studies should be made for predicting ducts 
more accurately so that one can rely on this phenomena. 
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APPENDIX: DERIVATION OF ELLIPSES FOR 

FRESNEL ZONES 



To derive the equations of the ellipsoid it is easier to change coordinates to x', y', z' 
by tilting the x coordinate by the angle of 8 in order to align the coordinate with R 
(see Figure 2.3). y' will be parallel to y and z' is obviously orthogonal to the x'y' plane 
[Ref. 3,4]. 

The expressions for the new coordinate systems x,' y', z' are [Ref. 4]: 



, ( r\ / U+t 2 \ • „ 

x = ( x- 1 c ° s p+ 1 2 - 2 lsin0, 


(A. 1 ) 


y' = y, 


(A.2) 


, ( r \ ■ a , ( Z \+ Z A a 

z — — ( x — ^ 1 sin0 + 1 2 — ^ 1 cos 8. 


(A.3) 


For given 6 the equation of the ellipsoid is: 




a 2 (y' 2 + z n ) + 6 2 x /2 = a 2 6 2 . 


(A.4) 


where 




r sec 8 + 6 

° 2 
1 


(A.5) 


b = -V6 2 + 2r6 sec 8. 

£ 


(A. 6) 


and 




tan«= Z! " Z ‘. 


(A.7) 



We are interested in the Fresnel zones in the xy plane transferring the ellipsoid Equa- 
tion A.4 back into x,y, z coordinates with z = 0 in the transformation equations, we 
get a family of ellipsoids: 



a 2 6 2 (b 2 cos 2 8 + a 2 sin 2 8 — c 2 ) / i \ 2 

^ 2 =(z--r-p) (6 2 cos 2 8 + a 2 sin 2 8 + a 2 y 2 ), (A.8) 

¥ cos^ 8 + a 1 sin 8 \ 2 / 
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where 



Z\ + Z2 

2 ’ 

c(a 2 — b 2 ) sin 9 cos 9 
b 2 cos 2 6 + a 2 sin 2 0 



(A.9) 

(A.10) 



Equation 2.5 represents an ellipse corresponding to a given 8 , the center of the 
ellipse is on x axis and is given by: 




1 - 



(^) 2 

(8/r 4- secf?) 2 — 1 



Vo = 0. 



(A.ll) 



By setting x = xq the semiminor axis is: 



Vi 



= ±b 




b 2 cos 2 6 + a 2 sin 2 9 ' 



(A.12) 



or 
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(2i + Zi ) 2 






[; + \A + er ) ; | 


2 

- 1 



By setting y = 0 the semimajor axis is derived: 



(A.13) 



x x = x 0 ± y t ===== . 

cos 2 # + l 2 sin 2 # 

Substituting Equations A. 5, A. 6, and A. 7 into Equation A. 14 yields: 



(A.14) 
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xo ± yj 
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(A. 15) 



For practical computations these equations can be simplified by 
Z \ , 22 and the path difference of 8 are much smaller than r. This yields 




assuming that 
to: 



(A.16) 
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